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On the geometry of broad emission region in quasars 
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ABSTRACT 

We study the geometry of the H/3 broad emission region by comparing the A/bh values 
derived from H/3 through the virial relation with those obtained from the host galaxy 
luminosity in a sample of 36 low redshift (z ~ 0.3) quasars. This comparison lets us 
infer the geometrical factor / needed to de-project the line-of-sight velocity component 
of the emitting gas. The wide range of / values we found, together with the strong 
dependence of / on the observed line width, suggests that a disc-like model for the 
broad line region is preferable to an isotropic model, both for radio loud and radio quiet 
quasars. We examined similar observations of the Civ line and found no correlation 
in the width of the two lines. Our results indicate that an inflated disc broad line 
region, in which the Carbon line is emitted in a flat disc while H/3 is produced in a 
geometrically thick region, can account for the observed differences in the width and 
shape of the two emission lines. 



Key words: 

lines 



galaxies: active - galaxies: nuclei - quasars: general - quasars: emission 



1 INTRODUCTION 

Super-massive black h oles (BHs) are found in virt ually 
all massive spheroids l|Kormendv fc Richstonel [l99^. In 
the Local Universe BH mass measurements can be per- 
formed th rough their imprint on the stellar kinematics (see 
iFerraresQ .2006, and references therein). The masses of the 
BHs are correla ted with some large s cale p r operties of their 



the underlying virial assumption (|Bottorf et al.l Il997^■ at 



least f or some broad em i ssion lines (e.g., Baskin fc Laoij 



12005^ ). iMcLure fc DunlopI l|200ll ) and iDunlop et all l|2003f) . 



host galaxies fr errarese fc Merrit t 
2000l: [Magorrian et al. 



;ale p i 
200d: 



11998 



Gebhardt et al,l 
Marconi fc Hunt 2003,; 



Graham fc Drive3l2007l '). The reader is referred to [Grahai 



2007 ) for an up-to-date review on this topic, and to iKin 



20051 ) and references therein for interpretative models of 



the scaling relations. 

In Type-1 AGNs the line emission from the gas inside 
the BH radius of influence is observed. If non-gravitational 
motions are neglected, the cloud velocity at a given ra- 
dius is fixed by the BH mass AIbh- Emission lines are 
Doppler broadened according to the gas motions. A sim- 
ple, isotropic mod el of the broad line reg ion (BLR) is usu- 
ally adop ted (e.g., Salviander et al. 120071 a nd the references 



the rein): [McLure fc DunlopI l|200a ). iDunlop erafl (|2003l ) 



Laor et al.1 (120061') found that a disc model is preferable, 
Decarli et al. I (2008') proved that inclination may explain 



and 
and 

some peculiar characteristics of the so-called Narrow Line 
Type-1 AGNs. On the other hand, some authors questioned 



and, more re cently, a number of reverberation mapping cam- 
paigns (e.g., Metzroth Onken fc Petersonll200^ : iBentz et al.l 
l2007l : ISergev et al.l I2007D found rough agreement between 
the virial estimate of the BH mass from H/3 width and 
the one based on the host galaxy luminosities, but the ob- 
served dispersions are significant. If the scatter is due to the 
assumed gas dynamical model, const r ains o n the BLR ge- 
ometry can be inferred. lOnken et aP l|2004 ) compared the 
virial estimates of Mbh in few, well studied nearby AGNs 
with the stellar velocity dispersion of their host galaxies. 
The large offset observed with respect to the relation ob- 
served in inactive galaxies suggested that an isotropic ge- 
ometry is not successful, but they could not put a better 
constrain on the gas dynamics due to the large uncertain- 
ties and the poor statistics. In radio loud quasars (RLQs), 
the width of the broad lines is roughl y correlated to the core- 
to-lobe power ratio index, Rc-i (e.g. IWills fc Brownelll986l : 
lBrothertonlll996l : IVestergaard Wilkes fc Barthelll2000D . This 
dependence is usually interpreted in terms of a flat BLR, 
given that -Rc-i is related to the inclination angle of the jet 
axis with respect to the line of sight. Whether a flat BLR 
model can be valid also for radio quiet quasars (RQQs) is 
still unclear. 



* roberto.decarli@mib.infn.it 



H/3 is the best studied emission line for low-redshift ob- 
jects, while MgllA2798 and C1VA1549 are often chosen for 
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Table 1. Sample properties and imaging data from literature. (1), (4)-(6): object name, redshift, radio loudness ('Q'=radio quiet; 
'L'=radio loud) and V apparent magnitude from Veron-Cetty & Veron (2006). (2)-(3): target coordinates from NED. (7): host galaxy 
absolute R magnitude corrected as described in section [3. II (8): estimated A/bh using equation [T] 



Object R.A. 
name (J2000.0) 
(1) (2) 



Dccl. 2 Radio 

{J2000.0) 

(3) (4) (5) 



V Ma log Mbh 

[mag] [mag] [Mq] 
(6) (7) (8) 
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higher redshift (e. g. McLure fc Jaryij 20021. Peterson et al.l 
|2004 iKaspi et al.l l2005l and 2007. iPeng et all l2006t ). since 
they fall in the optical range for z > 0.4 and z > 
1.6 r espectively. Mgll and H/3 widths are well correlated 
ISalvi andcr ct al. 2007,), while Civ line shows system- 
atic deviances from the H/3 values ijBaskin fc Laorl |2005| . 
IVestergaard fc Petersonll2006h . lEabita et al.l l|2006l ) studied 
Mbh derived from C iv width as a function of a Lbuigo- 
based Mbh for a sample of low-redshift quasars. They 
fo und a significantly better correlation than that reported 
bv lMcLure fc Dunlod (|200ll ) for H/3 data. 

In this paper we study the BLR geometry and gas dy- 
namics by comparing the Mbh values derived from H/3 broad 
emission with those obtained from the host galaxy luminos- 
ity in a sample of 36 low redshift {z ~ 0.3) quasars. The sam- 
ple is selected in order to provide similar numbers of RLQs 
and RQQs, so that conclusions on the geometry of the BLR 
can be drawn for bot h classes. Comparing our results with 
those of iLabita et al.l |200(|), based on Civ line, we sketch 
a picture of the dynamics of the gas around the BH. Note 
that since the dynamical model of the BLR is assumed to 
be independent of redshift, we also test the evolution of the 
BH-host galaxy relations. 



We define our sample in section [S] Data sources are 
summarized in section |3l Data analysis is described in sec- 
tion |4l We then discuss our results (section [5}. The sketch 
of the broad line region dynamics is presented in section |6l 
Throughout the paper we adopt a concordance cosmology 
with Ho = 72 km/s/Mpc, = 0.3 and = 0.7. We con- 
verted the results of other authors to this cosmology when 
adopting their relations and data. 



2 THE SAMPLE 



We selected all quasars in the IVeron-Cettv fc Veroiil (j2006l ) 
catalogue that have been imaged by i/S'r-WFPC2 (expo- 
sure time > 1000 s). The host galaxies are required to be 
elliptical, and therefore the bulge component practically co- 
incides with the whole galaxy. We considered all the quasars 
with z < 0.6, so that the H/3a486i-[0 ill] aa 4959,5007 lines 
are present in the optical spectra. For their observability 
from the Northern hemisphere, we selected only objects with 
5 > 0° . The entire sample thus consists of 53 targets. 12 of 
them have available spectra in the Sloan Digital Sky Survey 
(SDSS; see section [32J. Other 32 spectra were taken on pur- 
pose at the Asiago Observatory (see section [3^ . including 
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6 targets already observed by the SDSS. Two observed tar- 
gets were then excluded from our analysis: 0903+169 (RQQ) 
because of the low signal to noise ratio in the available 
spectrum, and 0923+201 (RLQ) because of the peculiar, 
composite profile of its broad lines, possibly due to the in- 
teraction with a ne arby galaxy (see iMcLeod fc Riek3ll994l : 
iBennert et al]|2007l ). Thus ~ 70 per cent of the sample was 
covered, including 16 RLQs and 20 RQQs. Table [1] summa- 
rizes the main properties of the objects in our sample. 



3 DATA SOURCES 

3.1 Host galaxy magnitudes and A/bh 



from the literature ( 


HooDer et al 


19971: Bovce et al. 


199S 


Kirhakos et al."l999|; 


PaEani et al. 


20031: iDunloD et all 


200S 



Labita ct al. 200g) o r converted from p ublished V or 
F702W luminosities (iBahcall et all 1 19971 : [Hamilton et all 



|2002| : iFlovd et al.l 20041) . Corre ctions for galactic extinction 
are from Schlegel et al.l (|l998l ). To perform colour and k- 



correction transformations we adopted an elliptical galaxy 
template (Mannucci et al. 2001), assuming that the host 
galaxies are dominated by old stellar population. The k- 
correction for an elliptical galaxy sd, z = 0.3 observed in the 
_R-band is 0.3 mag. The passive evolution of the host galax- 
ies foltows_^ressax|_jChiosi & Fagotto (1994), as discussed 
in iLabita et al.l l|2006l ) . Typical corrections for the passive 
evolution are ~ —0.3 mag. The resulting 7?-band absolute 

magnitudes of the host galaxies are given in t able U 

W e use the relationship obtained by iBettoni et al.l 
l|2003l l, corrected for the chosen cosmology, in order to esti- 
mate Mbh from the host galaxy luminosity: 



logA'/BH = -0.50 Mr - 2.60 



(1) 



where Mr is the absolute magnitude of the bulge component 
of the host galaxy. The rms of this fit is 0.39. Table [1] also 
lists the resulting Mbh values. 



3.2 Spectroscopic observations 

32 optical spectra were taken with the 1.82m Cima Ekar tele- 
scope at the Asiago Observatory. The Asiago Faint Object 
Spectrograph Camera was mounted in longslit spectroscopy 
configuration with grisms n. 4, 7 and 8, yielding spectral res- 
olutions of 7? ~ 300, 555 and 900 (2.10" slit) in the spectral 
range 3500-7800 A, 4300-6500 A and 6200-8050 A respec- 
tively (AA/pxl = 4.24, 2.10 and 1.78 A/pxl). At A ^ 5000 A 
the spectral instrumental resolutions are ~ 17, 9.1 and 5.5 
A, tight enough to distinguish prominent H/3 narrow emis- 
sion from the broad one. 

The standard IRAF procedure was adopted in the data 
reduction. The ccdred package was employed to perform 
bias subtraction, fiat field correction, image alignment and 
combination. Cosmic rays were eliminated by combining 3 or 
more exposures of the same objects, and applying crreject 
algorithm while averaging. When only one or two bidimen- 
sional spectra were available, we applied cosmicrays task 
in the crutils package. In order to prevent the task from 
altering the narrow component of emission lines, we masked 
the central region of our bidimensional spectra. The spectra 



Table 2. Sample objects spectroscopically observed at the Asi- 
ago Observatory. Redshifts and V magnitudes are taken from 
Veron-Cetty & Veron (2006). In the "Available spectra" column, 
Ax=Asiago + grism n.x; S=SDSS; U=HST-FOS archive data 
analyzed in Labita et al. (2006). Seeing is measured on corol- 
lary i?-band imaging. These images were not available for objects 
0100+020 and 3C48. 



Object 




Y 


Available 


Date 


Seeing 








spectra 




[arcscc] 




(2) 


(3) 


(4^ 


(5) 




0054-1-144 


0.171 


15.70 


A4,S 


13/09/06 


1.4 


0100-1-020 


0.393 


16.39 


A4,H 


18/09/07 


n/a 


0110-1-297 


0.363 


17.00 


A4 


09/12/05 


2.1 


0133-1-207 


0.425 


18.10 


A4,H 


27/11/06 


2.0 


3048 


0.367 


16.20 


A4 


11/12/05 


n/a 


0204-1-2916 


0.109 


16.80 


A4 


15/10/06 


1.7 


0244-1-194 


0.176 


16.66 


A4 


16/10/06 


1.7 


0624-1-6907 


0.370 


14.20 


A4,H 


27/10/06 


2.6 


07546-1-3928 


0.096 


14.36 


A4 


24/12/05 


1.9 


US1867 


0.513 


16.40 


A4,S,H 


27/11/06 


1.3 


0944.1-1-1333 


0.131 


16.05 


A4 


25/04/06 


2.4 


0953-1-415 


0.234 


15.30 


A4 


07/02/06 


3.1 


1001-f291 


0.330 


15.50 


A4 


19/02/06 


2.4 


1004-1-130 


0.240 


15.20 


A4,S 


07/03/06 


2.2 


1058-f 110 


0.423 


17.10 


A4,S 


26/04/06 


2.0 


1100-1-772 


0.315 


15.72 


A4 


08/03/06 


2.7 


1202-1-281 


0.165 


15.60 


A4,H 


20/02/06 


2.3 


1216-1-069 


0.331 


15.65 


A8,H 


12/04/07 


1.9 


Mrk0205 


0.071 


15.24 


A7,H 


11/04/07 


1.6 


1307-1-085 


0.155 


15.10 


A7,H 


24/04/07 


1.6 


1309-1-355 


0.184 


15.64 


A4,H 


20/02/06 


2.4 


1402-1-436 


0.323 


15.62 


A4,S 


07/03/06 


3.6 


1425-1-267 


0.366 


15.68 


A4,H 


20/02/06 


2.9 


30323.1 


0.266 


16.70 


A7,H 


24/04/07 


1.4 


1549-1-203 


0.250 


16.40 


A4 


23/04/06 


2.0 


1635-H119 


0.146 


16.50 


A4 


31/05/06 


2.8 


1821-1-643 


0.297 


14.10 


A4,H 


15/12/06 


2.1 


2141-1-175 


0.213 


15.73 


A4,H 


15/12/06 


1.9 


2201-1-315 


0.295 


15.58 


A4,H 


15/12/06 


1.5 


2247-1-140 


0.235 


16.93 


A4,S,H 


12/09/06 


1.5 



extraction, the background subtraction and the calibrations 
both in wavelength and in fiux were performed with doslit 
task in kpnoslit package, using a Hg-Cd and Ne lamps and 
spectrophotometric standard stars as reference. Wavelength 
calibration residuals are around 0.17 A (sub-pixel), thus im- 
plying a negligible (< 1 per cent) error on redshift estimates. 
Absolute calibration of spectra was corrected through the 
photometry of field stars, by comparing corollary imaging 
with Johnson's R and V filters to the magnitudes published 
in the U.S. Naval Observatory catalogue. The uncertainty 
in the flux calibration is 0.1 mag. Galactic extinction was 
accounted for according to Schlegel, et al. (1998) , assuming 
Rv = 3.1. 

Table [2] summarizes the observed targets. 

The SDSS Data Release 5 (Adelman-McCarthy, et al., 
2007) provides spectra for 12 quasars in our sample. SDSS 
spectra have R ~ 2000 and a spectral range between 3800 
and 9000 A. Uncertainties on wavelength calibration amount 
to 0.05 A, while flux calibration formal errors account to 5 
per cent. We re-observed six of these objects (0054+144, 
US1867, 1004+130, 1058+110, 1402+436, 2247+140) in or- 
der to perform a comparison. Both the speciflc fluxes at 5100 
A and the H/3 broad line widths are in good agreement. Due 
to their better spectral resolution, we will consider only the 
SDSS spectra of these objects in our analysis. 
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4 THE VIRIAL DETERMINATION OF Mbh 

In the virial assumption, if the velocity u of a particle orbit- 
ing at a certain radius R around the BH is known, the mass 
of the BH is simply: 



2 = 0.297 



Mb 



G 



(2) 



In Type-1 AGNs expression ([2]) can be evaluated at the char- 
acteristic radius of the BLR, R = -Rblr, which can be mea- 
sured almost directly with the re v erberation mappin g; tech- 
nique (jBlandford fc McKedlTosj ) . iKaspi erall l|2000l ) found 
that the nuclear monochromatic luminosity, \L\, is corre- 
lated with -Rblr by \Lx oc i?BLHi with 7 depending on the 
considered emission line. The broad line width is used to 
estimate «blr = «(-Rblr)- According to the adopted dy- 
namical model, 

^^BLR = / • FWHM (3) 

where / is a geometrical factoiQ and the FWHM is expressed 
in velocity units. If the gas moves isotropically, a gaussian 
line shape is observed, with / ~ If a rotational com- 

ponent is present, an axial symmetry is introduced, and / 
depends on the inclination angle •&: 



(4) 



/ = I 2ci sini9 + "^^2 



where •& is the angle between the line of sight and the ro- 
tation axis, and ci and C2 are parameters accounting for 
the importance of the disc and isotropic components resgec- 
tively. The r e ader is referred to (McLure fc Dunlopl ([200 If). 
iLabita et all l|2006l ). ICoUin et al.l (120061 ) and iDecarh 
( 2008h for detailed discussions on this topic. Since / value 
is generally unknown, one defines the Virial Product as: 



VP = AfBH ■ = G-^RBLR{VWmAf 



(5) 



corresponding to Mbh if / is taken unitary. Thus, an esti- 
mate of / requires measures of the monochromatic luminos- 
ity \L\, of the line width and independently of Mbh- 



4.1 The monochromatic luminosity 

The continuum luminosity was calculated as follows: Rest- 
frame spectral continua were fitted with a power-law (see 
figure [T]) . The key-point is to avoid contaminations due to 
various emission features, especially to broad Fell bands. 
Typical rest-frame fitted regions are around 2610, 3030, 
4030, 4190, 5080 A and in the 5450-5720 A range. The ~ 
5100 A region is generally free from contaminations and can 
be easily fitted. The resulting fit function was computed at 
5100 A, thus providing AFa(5100 A). We corrected for the 
host galaxy contamination by computing the fraction of its 
flux within the slit, on the basis of the nucleus and host 
galaxy luminosities and the effective radius estimates avail- 
able from the literature , and following the recipe adopted by 
ISbarufattT et all (|2006l '). The average correction is 0.06 dex. 



^ The reader should noti ce that different defini t ions for / 
are available in li terature: IVestergaard fc Petersonl ( l2006tl and 
ICoIIin et al.l (I2OO6I ') use Mbh = / ■ G-1Rblr(FWHM)2. 




Figure 1. The rest-frame spectrum of 1821-1-643. The 
upper panel shows the observed spectrum with the fit- 
ted power-law for the continuum. The lower panel shows 
the residual after continuum subtraction. Main emission 
lines and the atmospheric absorption features are also la- 
beled. All the spectra are available in electronic format at 
www. dfm. uninsubria. it / astro / caqos / index, html . 



4.2 Line width measurements 

The observed H/3 broad component is usually contaminated 
by other spectral features, in particular by blended Fe 11 mul- 
tiplets, the H/3 narrow component and the [O ill] AA4959, 5007 
lines. 

Strong Fe 11 emissions are commonly detected in quasar 
spectra at 4400-4750 and 5150-5450 A, and weaker blended 
features are found at 4800-5000 A. A common practice 
to remove this contamination (e.g., Boroson fc Green 19921: 
Marziani erahl l2003l : ISalviander et all l2007l : iMcGill et all 
20071 ) is to adopt the spectrum of I ZwOOl as a template 
of the Fell emission, due to the intensity and narrowness 
of the Fell lines. Since the relative intensities of the vari- 
ous Fen features may differ from a q uasar to another (see, 
for example, table 7 in |Phillips|[l97i ). we preferred a more 
conservative approach: we modeled the Fell emission as a 
simple power-law fitted at ~ 4750 and ~ 5100 A, and sub- 
tracted it from the observed spectrum. The reliability of this 
procedure is discussed in appendix lAll 

Since most of our objects have the H/3 red wing contam- 
inated by [Fevill]A4S94, FellAA4924,50i8 and [O ill] aa4959, 5007 
lines, a reliable study of the broad line asymmetries is ex- 
tremely hard to achieve, and is strongly dependent on the 
procedure adopted in removing these contaminations. We 
preferred to set the same peak wavelength for both the gaus- 
sian curves, thus neglecting line asymmetries. In appendix 
A2l we discuss how the use of a different fitting function, sen- 
sitive to the asymmetries in the line profile, does not affect 
the estimates of the line width in a significant way. 
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Table 3. Continuum luminosity parameters. (1): object name. 
(2): object measured redshift. (3): data set. 'A' refers to Asiago 
data, the number referring to the adopted grism; 'S' refers to 
SDSS data. (4): log i^>,{5100A), computed as described in section 
14.11 (5): final logAL;^ [erg/s]. (6): computed log/?BLR, [cm]. 



Object 


z 


Data 


log Fa 


log ALa 


log i?BLR 


name 




set 


[crg/s/cm'^ / A] 


[crg/s] 


[cm] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


0054+144 


0.171 


g 


-14.9 


44.7 


17.2 


0100+0205 


0.393 


A4 


-15.6 


44.9 


17.3 


0110+297 


0.363 


A4 


-15.5 


44.7 


17.2 




0.424 


A4 


-15.6 


44.7 


17.2 


3C48 


0.369 


A4 


-14.9 


45.4 


17.6 


0204+292 


0.110 


A4 


-15.1 


44.2 


16.8 


0244+194 


0. 174 


A4 


-15.1 


44.7 


17.2 


nf\9d-i-f^QD7 


0.370 


A4 


-14.5 


46. 1 


18.1 




0.096 


A4 


-14.2 


44.7 


17.2 


TTQl ftfi7 


515 


g 


-15 3 


45 3 


17 6 


±-1" lOOO 


134 


A4 


-14 8 


44 5 


17 1 




0.235 


A4 


-14.8 


45.6 


17.8 




330 


A4 


-14 9 


45 1 


17 5 


1 nn4-i-i '^0 


0.241 


s 


-14.7 


45.2 


17.5 




0.423 


g 


-15.9 


44.6 


17. 1 


1100+772 


0.311 


A4 


-14.4 


45.8 


17.9 


1150+497 


0.334 


S 


-15.6 


44.6 


17.1 


1202+281 


0.165 


A4 


-15.3 


44.6 


17.1 


1216p069 


0.331 


A8 


-14.7 


45.6 


17.8 


Mrk0205 


0.071 


A7 


-14.7 


44.5 


17.0 


1222+125 


0.412 


S 


-15.6 


44.8 


17.2 


1230+097 


0.416 


S 


-15.2 


45.2 


17.5 


1307p085 


0.155 


A7 


-15.0 


45.3 


17.6 


1309+355 


0.184 


A4 


-14.7 


44.4 


17.0 


1402+436 


0.323 


S 


-14.7 


45.4 


17.7 


1425+267 


0.366 


A4 


-15.4 


45.2 


17.5 


1444+407 


0.267 


S 


-14.9 


45.1 


17.4 


1512+37 


0.371 


S 


-15.3 


45.0 


17.4 


3C323.1 


0.266 


A7 


-15.2 


45.2 


17.5 


1549+203 


0.253 


A4 


-15.2 


44.8 


17.3 


1635+119 


0.148 


A4 


-15.4 


42.2 


15.5 


3C351 


0.372 


S 


-14.7 


45.6 


17.8 


1821+643 


0.297 


A4 


-14.2 


46.1 


18.1 


2141+175 


0.211 


A4 


-15.0 


45.0 


17.4 


2201+315 


0.295 


A4 


-14.9 


45.4 


17.7 


2247+140 


0.235 


S 


-15.3 


44.6 


17.1 



The fit procedure was preferred to the width measure- 
ment directly a pplied to the observ ed data (without any fit; 
see for example ICollin et al.ll200^ since: 1) it is applicable 
also to low signal-to-noise spectra; 2) it does not require 
an accurate modeling of the narrow component; 3) it is reli- 
able even in the largest tails, where contaminations by other 
emission or absorption features may be relevant. 

We derived the FWHM and the second moment of 
the line, aune, from the fitted profile. Both FWHM and 
criine are corrected for instrumental spectral resolution. 
The ratio between FWHM and aunc is used to study the 
shape of the line, as discussed in section 15.11 According to 
ICoUin et aTl (|2006l ). criine could be preferred to FWHM as a 
line width indicator, since it is strongly dependent on the 
line wings, i.e., to the kinematics of the innermost clouds. 
On the other hand, aunc is very sensitive to contamina- 
tions and to the adopted Fell subtraction. Therefore, we 
will use the FWHM in the Mbh estimates, and consider 
(Tiine only in the study of the line shape. Figure [2] offers 
an example of the fitting procedure, applied to 1821-1-643. 
Plots for the whole sample are available in electronic form 
at www.dfm.uninsubria.it/astro/caqos/index.html. Our 
H/3 width estimates are listed in table [l] Typical uncertain- 




4700 4800 4900 5000 5100 

m 

Figure 2. The rest-frame spectrum of 1821+643 in the 4660-5160 
A range. The histogram is the observed spectrum, after contin- 
uum and Fell subtraction, as described in the text. The solid, 
smooth line is the fitted curve. The residual after fit subtraction 
is plotted with a dashed line. A vertical, dotted line shows the 
expected H/3 wavelength (4861 A). 



ties in the FWHM values due to the fit procedure are ~ 10 
per cent of the line width. 



5 DISCUSSION 

5.1 H/3 line width and shape 

The mean value and rms of the H/? FWHM distribution are: 

< FWHM (H/3) >Aii= 5050 ± 2170 km/s, 

< FWHM (H/3) >RLQs= 6100 ± 2110 km/s, 

< FWHM (H/3) >RQQs= 4210 ± 1870 km/s. 

A different distribution of FWHM is observed between RLQs 
and RQQs, in the sense that the former ones show, on aver- 
age, wider lines than the latter. This difference may be in- 
trinsic, the velocity of the gas in RLQs being actually larger 
than in RQQs, independently on the BLR geometry. Indeed, 
as we will notice later, a bias towards the high-mass end of 
the Mbu distribution occurs for RLQs (partially related to 
the Malmquist bias, since the average redshift of our sam- 
ple RLQs is higher than that of RQQs). Otherwise, in the 
prospective of a disc-like broad line region, different average 
inclination angles may account for the different distributions 
in the FWHM values. This may be achieved assuming that, 
while RLQs have 0° < < 50°, the RQQs are biased to- 
wards lower inclination angles, e.g. 0° < i9 < 40°. Such a 
bias has been alrea dy hypothesized by some authors (e.g., 
iFrancis et al 



for example. 



20001). but its occurre nce is still debated (see. 



Kotilainen et al.ll2007l ') 
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Table 4. The estimated FWHM and an^^ values for H/3. Asiago 
spectra are labeled with 'Ax' in column (2), where x is the grism 
number. SDSS data are labeled with 'S'. 



Object 


Data 


FWHM 




FWHM 


^linc 


name 


set 


[A] 


[A] 


[km/s] 


[km/s] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


0054+144 


S 


134 


57 


8220 


3520 


0100+0205 


A4 


79 


35 


4870 


2160 


0110+297 


A4 


100 


40 


6150 


2440 


0133+207 


A4 


143 


53 


8850 


3260 


3C48 


A4 


65 


39 


4010 


2370 


0204+292 


A4 


145 


67 


8960 


4110 


0244+194 


A4 


76 


36 


4680 


2220 


0624+6907 


A4 


59 


32 


3630 


1940 


07546+3928 


A4 


50 


28 


3110 


1750 


US1867 


S 


37 


31 


2280 


1900 


0944.1+1333 


A4 


52 


30 


3230 


1850 


0953+415 


A4 


58 


47 


3610 


2890 


1001+291 


A4 


32 


26 


1980 


1610 


1004+130 


S 


98 


50 


6010 


3060 


1058+110 


S 


121 


48 


7460 


2950 


1100+772 


A4 


113 


58 


6980 


3550 


1150+497 


S 


63 


28 


3870 


1750 


1202+281 


A4 


84 


57 


5190 


3510 


1216+069 


A8 


83 


50 


5110 


3080 


Mrk0205 


A7 


53 


34 


3270 


2100 


1222+125 


S 


122 


56 


7530 


3430 


1230+097 


S 


77 


41 


4770 


2540 


1307+085 


A7 


57 


27 


3520 


1670 


1309+355 


A4 


72 


44 


4430 


2720 


1402+436 


S 


45 


33 


2760 


2030 


1425+267 


A4 


131 


73 


8090 


4530 


1444+407 


S 


44 


27 


2690 


1650 


1512+37 


S 


144 


52 


8910 


3210 


3C323.1 


A7 


77 


37 


4760 


2290 


1549+203 


A4 


30 


20 


1880 


1230 


1635+119 


A4 


92 


46 


5700 


2870 


3C351 


S 


150 


63 


9260 


3870 


1821+643 


A4 


78 


34 


4820 


2090 


2141+175 


A4 


84 


36 


5180 


2230 


2201+315 


A4 


49 


28 


3020 


1720 


2247+140 


S 


52 


22 


3180 


1330 



The comparison between H/3 FWHM and crunc illus- 
trates some general properties of the line shape (see figure 
[3|. A single gaussian has FWHM/(Jiine= V81n2 ^ 2.35. The 
single gaussian case (upper dashed line) represents an up- 
per limit for FWH M. Only few H/? da ta have FWHM<(Jiine 
(lower dashed linel. lCollin et al.l l|2006h suggested a bimodal- 
ity in the FWHM vs (June relation, with a break when aiine 
> 2000 km/s. We argue that such a behaviour is mainly 
due to the fit procedure adopted by those authors: When 
the H/3 width is estimated directly on the observed spec- 
trum, (Tiino cannot be integrated up to infinity, because of 
He IIA4686, [Ar iv]aa47ii,4740 and Fell contaminations. Inte- 
gral truncation leads to underestimates of the largest line 
widths: For a single gaussian curve, the deviation is sig- 
nificant when (Tiine^ 0.5 times the width of the truncation 
interval. Typically, H/3 can be studied only in the first 80 A 
bluewards. That means, crunc > 2500 A are underestimated. 
In disagreement with Collin's results, our fit-based FWHM 
to cTiine ratio is found to be constant all over the observed 
values of aunc- No systematic difference is reported in the 
FWHM/CTiine ratio of RLQs and RQQs. 
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Figure 3. The FWHM as a function of a^^^ for H/3 (filled sym- 
bols) and C IV (empty symbols). Circles refer to the RLQs, squares 
to the RQQs. The dashed lines refer to FWHM= 2.35 tTuno and 
FWHM=(Tii„c. A typical error box is also plotted. H/3 and Civ 
data clearly fill different regions of the plot. No significant differ- 
ence in the FWHM/criine ratio is observed as a function of the 
radio loudness. 



5.2 Broad Line Region radii 

Kaspi et al. (2000) found that the radius of the broad 
line region, as estimated by the reverberation mapping, is 
related to the continuum monochromatic luminosity AL^- 
An increasing bo dy of measurements is now available for 



H/3 time lags (e.g . iKaspi et aL 2000|. Suganuma et al.ll2006l , 
iBentz et al.ll2006l ) . Following iKaspi et all l|2005h we adopt : 



°< X \ 0-87+0.07 

,„.„ / ALa(5100^) \ 
10 light -days ^—"-024) ■ I ig** erg/s J 



i?SLfl (H/3) 



(2.OO: 



(6) 



The average characteristic radius for H/3 is: 

< logi?BLR (H/3) [cm] >Aii= 17.35 ± 0.44 

< logiiBLR (H/3) [cm] >RLQs= 17.38 ±0.28 

< logi?BLR (H/3) [cm] >RQQs= 17.33 ±0.55 
where the error is the standard deviation. 



5.3 Redshift dependence of the Afen— ibuige 
relation 

Our work is centered on the comparison between the Virial 
Products (equationj5|_and__the Men evalu ated from t h e host 
galaxy luminositv. IWoo et al.l (|2006t ) and lTreu eraP (|2007[ ) 
proposed that the Men-bulge relations change significantly 
even at z ~ 0.36. On the other hand, iLauer et al.1 (120071 ) 
showed that such a result is probably due to a statistical 
bias. Owing to the steepness of the bright end of galaxy lu- 
minosity function, very high mass black holes are more com- 
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monly upper residuals of the MeH-host luminosity, rather 
than being hosted by correspondingly very massive galax- 
ies. Thus, when a luminosity cutoff is adopted (typically 
due to sensitivity limits, when observing at high redshift), 
the Mbh expected from MBH-ibuigc relation tends to be 
lower than the real one. Since the geometrical factor is sup- 
posed to be independent of redshift, we can directly test 
the evolution of the Bettoni relation (equation [ij checking 
the redshift dependence of VPs, host absolute magnitudes 
and their ratios. We consider here only the objects the host 
galaxies of which have similar luminosities, namely, the (15) 
objects with — 23 > Mr > — 24 mag, since they are roughly 
well distributed along the considered redshift range. While 
an overall slight increase in the VPs is found with redshift 
(possibly due to Malmquis t bias), data dispe rsion largely ex- 
ceeds the effect reported in IWoo et al\\200 &) and Tr eu et al.l 
our objects being consistent with a no-cvolution sce- 
nario. A low-er scat ter is observed for C iv data taken from 
iLabita et al.l (|2006l ) (see section 15. 5p . Applying the same ar- 
gument, no significant redshift dependence is found in the 
VP-to-host galaxy luminosity ratios, the probability of null 
correlation exceeding 30 per cent. 



5.4 Spectroscopic Virial Products vs Imaging 
A/bh estimates 

We are now ready to compare the Virial Products to the BH 
masses as estimated from the host galaxy luminosity. VPs 
and / values are listed in table [S] Typical uncertainties on 
VPs and Mbh are ~ 0.4 dex, due mainly to the scatter in 
the -Rblr-ALa and MBH-i/buigo relations. The comparison is 
shown in figure|4] Our data show no correlation: The proba- 
bility of non-correlation is ~ 40 per cent, with a Spearman's 
rank coefficient of 0.20 and a residual standard deviation of 
~ 0.50 dex. The mean value of / is 1.6± 1.1. The dispersion 
in our data reflects the values obtained in the same way by 
iDunlop et all l|2003l ) (see their figure 13). 

We study now the geometrical factor of our sample 
quasars separately according to radio loudness. Indeed, in- 
dications that RLQs may have flat BLRs have been already 
reported (in partic ular, the rough dependence of Rc-i on 
the FWHM of H/3:rwills fc Browne|ll986l : lBrothertonlll996l : 
IVestergaard Wilkes fc Barthelll200ol ). while not so much is 
known about RQQs. In RLQs, the / factor is found to be 
strongly dependent on the FWHM of the lines, as shown in 
figure [SI upper panel. We note that such a relation show 
significantly less dispersion than the FWHM -_Rc-i pub- 
lished for similar samples (e.g., see figure 4 in iBrothertonI 
1 19961 ). Some authors proposed the existence of two differ- 
ent po£ulations_wM^ and average 
/ (|Sulentic et al.l |200G| : ICollin et al1l2006l ). Our data con- 
firm this result, even if a continuous trend rather than a 
strict bimodality seems to be present. All the objects with 
FWHM< 5000 km/s have / > 1, the value of / rapidly de- 
creasing with FWHM. According to equation|31 such a trend 
would be expected only if / is not fixed, given wblr- This 
reinforces the idea that RLQs have disc-like BLRs. 

A similar picture is observed for RQ Qs, even if a larger 
dispe rsion is found (see, for comparison, iMcLure fc Dunlog 
|2002| ). The occurrence of two populations is still clear: at 
FWHM< 4000 km/s all but two targets have / > 1, while 



Table 5. Virial Products and geometrical factors for optical data. 
(1): object name. (2): object measured redshift. (3): Virial Prod- 
ucts derived from H/3 width and 5100 A monochromatic lumi- 
nosity. Error bars are mainly due to the intrinsic scatter in the 
^BLR^-^^A relation (~ 0.4 dex). (4): / values. 



Object 


z 


log VP (U/3) 


/ 


name 








(1) 


(2) 


(3) 


(4) 


0054+144 


0.171 


8.9 


0.7 


0100+0205 


0.393 


8.6 


0.8 


ni 1 n 1 ori*7 
OiiU+^y ( 


0.363 


8.6 


1 . 1 


Oioo+^U f 


0.424 


9.0 


0.6 


3C48 


0.368 


8.7 


3.0 


0204+292 


0.110 


8.6 


1.1 


AO A A 1 1 A /I 

Uii44+iy4 


0.174 


8.4 


1.1 


0624+6907 


0.370 


9.1 


1.8 


07546+3928 


0.096 


8.0 


4.6 


US1867 


0.515 


8.2 


2.6 


0944.1+1333 


0.134 


7.9 


3.4 


0953+415 


0.235 


8.8 


0.7 


iooi+zyi 


0.330 


7.8 


4.4 


1 A A A 1 1 O A 


0.241 


9.0 


0.9 


1 A C O 1 1 1 A 

iUoo + iiU 


0.423 


8.7 


0.8 


1100+772 


0.311 


9.5 


0.6 


1150+497 


0.334 


8.1 


3.1 


1202+281 


0.165 


8.4 


1.2 


1216+069 


0.331 


9.1 


0.5 


Mrk0205 


0.071 


8.0 


2.2 


1222 + 125 


0.412 


8.8 


1.0 


1230+097 


0.416 


8.8 


1.7 


1307+085 


0.155 


8.5 


0.7 


1309+355 


0.184 


8.1 


2.2 


1402+436 


0.323 


8.5 


1.5 


1425+267 


0.366 


9.2 


0.7 


1444+407 


0.267 


8.2 


1.8 


1512+37 


0.371 


9.2 


0.7 


3C323.1 


0.266 


8.8 


1.1 


1549+203 


0.253 


7.5 


2.4 


1635+119 


0.148 


6.9 


0.6 


3C351 


0.372 


9.6 


0.5 


1821+643 


0.297 


9.4 


1.3 


2141+175 


0.211 


8.7 


1.3 


2201+315 


0.295 


8.4 


3.2 


2247+140 


0.235 


8.0 


2.7 



only few targets with FWHM> 4000 km/s have / > 1. This 
rules out that the BLR is isotropic even in RQQs. 

We simulated the /-FWHM relation in the hypothe- 
sis of a thin disc BLR (figure [SJ. We assumed a gaussian 
distribution of the uncertainties for FWHM, log VP and 
log Mbh , mimicking the uncertainties in the adopted fitting 
techniques and scaling relations. We then assumed that all 
the quasars have purely disc-like BLRs, with a fixed rota- 
tional velocity wblr- Fitting our data with a hyperbole (see 
eq. [3]), we found i;blr ~ 8000 and 6000 km/s for RLQs and 
RQQs respectively. The variable was let free to vary from 
to i9max. The angle i?max was fixed to 40° and 50° for RLQs 
and RQQs respectively, in order to match the median value 
of the / observed distributions. The simulated values are 
overplotted to the measures presented in figure [5] for a com- 
parison. We also plotted the two lines corresponding to the 
expected values of / if uncertainties were negligible. Both 
RLQs and RQQs are well described with a disc model of the 
BLR, with the former ones showing larger vblb. than the 
latter ones. This difference cannot be explained in terms of 
a different range of i?, and - as we noted in section 15.11 - 
may be the effect of a selection bias, in the sense that the 
RLQs in our sample have higher average Mbh than RQQs. 
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Figure 4. Comparison between spectroscopic Virial Products 
based on H/3 and C IV and Mbh estimates based on the bulge 
luminosity. Symbols are the same as in figure |3] The dashed line 
is the expected VP as a function of Mbh for a-n isotropic model. 
The values expected for a geometrically thin disc model with dif- 
ferent inclination angles (labeled in degrees) are plotted as dotted 
lines. 



Figure 5. Comparison between the geometrical factor, /, and 
the FWHM for RLQs (upper panel) and RQQs (lower panel). 
Symbols are the same as in figure |3] The isotropic case is shown 
as a dashed line, while the values of / expected for a thin disc 
model with diff'erent inclination angles (labeled in degrees) are 
plotted as dotted lines. 



Concerning the different FWHM distributions of RLQs and 
RQQs, our simulation cannot rule out a dependence on the 
adopted range of i9, as suggested by the the estimates of 
i^max based on the distributions of /. Extending this tech- 
nique on a larger sample is needed to properly address this 
topic. 

5.5 Comparison between Civ and H/3 lines 

We now want to compare the H/9 properties wit h those of the 
Civa 1549 line. We will refer to the C iv data in lLabita et al.l 
who used the C iv line in order to measure the VPs 
of 29 low-redshift quasars observed with HST-FOS. 16 out 
of 29 objects in that work are also present in our sample. 
The average C iv FWHM and its standard deviation are: 

< FWHM (CIV) >Aii Labita sample 

= 4200 ± 1200 km/s, 

< FWHM (CIV) >common sample^ 4030 ± 1200 km/s. 

C IV FWHMs show smaller dispersion and a smaller mean 
value than H/3 ones. This is remarkable: the C iv line re- 
quires higher ionization potential than H/3. If the virial hy- 
pothesis is valid and a simple photoionization model is as- 
sumed, C IV emission should show lower radii and, accord- 
ingly to equation [21 higher velocities. In section |B] we will 
show how this point can be interpreted in terms of differ- 
ent thickness of the BLR disc at C iv and H/3 radii. Figure 
\7\ shows the comparison between C iv and H/3 FWHMs for 
the objects common to both our and Labita's studies. No 
correlation is apparent (the probability of non-correlation 




2000 4000 6000 8000 10000 12000 

FWHM [km/s] 

Figure 6. The /-FWHM relation for RLQs (upper panel) and 
RQQs (lower panel), as described in the text. Dots refer to the 
simulated values. The solid lines show the expected / assuming 
a thin-disc geometry with dblr = 8000 and 6000 km/s for RLQs 
and RQQs respectively. 
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Figure 7. Comparison between H/3 and C IV estimates of FWHM. 
Circles mark RLQs, squares refer to RQQs. A typical error box 
is also plotted. The dashed line is the one-to-one relation. No 
correlation is observed. 



being up to ~ 4 5 per cent), as alread y noticed for dif- 
ferent samples by Baskin fc Laorl (|2005l ) (81 quasars) and 
IVestergaard fc Peterson (200^ ) (32 quasars; see their figure 
10). We argue that the FWHM values of the two lines are 
intrinsically different. 

As in the H/3 case (figure [3|, we now consider the com- 
parison between C iv FWHM and ffune- Even if we sampled 
only a small range of aunc values, the data from the two 
lines clearly fill different regions of the (FWHM, aunc) plane. 
Since the fitting procedure is similar, we can rule out a sys- 
tematic effect due to the width estimate algorithm. Thus, 
we argue that C iv and H/3 broad lines have intrinsically dif- 
ferent shapes, the C iv departing more from the isotropic 
(gaussian) case than H/3. 

The continuum monochromatic luminosity at 1350 A is 
also considered for the objects belonging to the sample of 
iLabita et al.l (j2006l ) . We use this information in order to esti- 
mate -Rblr for C IV d ata, by means of th e radius-luminosity 
relation published bv lKaspi etahl (|2007l ): 



i^BLR(CIV) 

10 light — days 



(0.24±0.06)- 



XL 



.(1350A) y 



10-^3 erg/s J 



(7) 



This is the first available relation based on reverberation 
mapping studies of the C iv lines. However, it is based only 
on 8 sources, thus the slope and offset of the relation shall 
need some tuning when other data will be available. The 
average -Rblr(Civ) result: 

< logi^BLR (CIV) [cm] >A11 Labita sample 

= 17.19 ±0.31 

< log^BLR (CIV) [cm] >c ommon sample = 17.19 ±0.27 

The -Rblr(Civ) is found to be systematically (~ 1.7 
times) smaller than i?BLR(H/3), but we warn the reader that 



the scatter due to the -Rblr - ALa relations is severe. Val- 
ues of -Rblr(Civ) smaller than -Rblr(H/3) are consistent 
with simple photoionization models: The C iv line is a higher 
ionization line than H/3, thus it should be emitted in an in- 
ner region. Such a difference, if confirmed, should be taken 
into account when referring to previous works where the 
J?blr(H/3) was used as a surrogate for 7?blr(Civ). 

Following the steps traced in section 15.41 we estimate 
the VPs also for the C iv data. Carbon VPs are well corre- 
lated with imaging BH mass estimates (see figure 3)) : The 
probability of non-correlation is < 0.1 per cent with a Spear- 
man's rank coefficient of 0.63 and a residual standard devi- 
ation of ~ 0.33 dex, comparable with the dispersion in the 
luminosity-radius relations. Even considering only the com- 
mon sample, the probability of non-correlation for C iv and 
H/3 data is < 0.1 and ~ 60 per cent respectively. The mean 
value of / for Civ data amounts to 2.40 ±0.16. This value is 
1.75 t imes larger than the value obtained by iLabita et all 
l|2006l ) on the same data, since they adopted a different 
^BLR - ALa relation (|Pian Falomo fc Treves! [20051 ) , which 
provided an estimate of the H/3 broad-line radius rather than 
the C IV one. The / dependence on FWHM is observed also 
in Civ data (see figures [S] and O . No / value deduced from 
C IV line is found to be consistent with an isotropic model 
of the BLR. The /-FWHM plot for Civ-based data is fully 
consistent with a disc-like BLR, both for RLQs and RQQs. 
In the average value of /, the C iv line appears farther from 
the isotropic case than the H/3, both being inconsistent with 
isotropy. 

The formally better correlation of imaging-based BH 
mass estimates with VPs for C iv line rather than 
for H/3 suggests that the former could be preferable 
as a mass indicator than H /3, as already proposed by 
IVestergaa rd fc PetersonI l|2006l ) and referenc es therein. If 
confirme d, this disagrees with recent claims (B askin fc Laorl 
l2005i : ISulentic et al.ll2007[ ). in which the average blueshift of 
the high-ionization lines with respect to the rest frame of the 
host galaxy was inte rpreted in terms of gas outflow (but see 
[Richards et aLir2002l for a different interpretation). A larger 
set of quasars with independent estimates of both VPs and 
Mbh, sampling a wider parameter space, is required in order 
to better address this topic. 



6 A SKETCH OF THE BLR DYNAMICS 

We discuss how three simple models of the BLR dynamics 
can account for or contrast with the observed H/3 and C iv 
line shapes and widths. 

Isotropic model: 
Up to now, since the geometry of the BLR is poorly un- 
derstood, an isotropic model has been commonly adopted 
as a reference. As mentioned in section 15.11 if the BLR is 
dominated by isotropic motions, with a Maxwellian veloc- 
ity distribution, the geometrical factor is / = V3/2 and 
the FWHM/criine ratio is %/8ln2 ^ 2.35. We find that H/3 
FWHM/(Tiinc ratios are closer to the isotropic case than C iv 
ones. All the geometrical factors derived from C iv and most 
of those from H/3 exceed unity, in disagreement with the ex- 
pected value for this model. Moreover, the isotropic model 
does not explain the observed /-FWHM relation. 
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Geometrically thin disc model: 

If a disc model is adopted, / depends on 3 free parameters, 
namely i?, ci and C2 in equation |4] Here we assume that 
the disc is geometrically thin, that is, C2 tends to zero, ci 
tends to 1 and / ~ (2sm'd)~^ . In the ligh t of the unified 
scheme of AGNs l|Antonucci fc Miller! llOSa ), is supposed 
to vary between i?min, fixed by the presence of a jet (if any) 
and i9max, given by the angular dimension of the obscuring 
torus. Even if these angles are poorly constrained, reason- 
able values are 0° < i^min < 10° and 30° < iJmax < 55° 
fo r Type-1 AGNs (se e for e xample the statistic al approach 
in iLabita et al.ll2006l and in iDecarh et all 120081 ). Since the 
line profile depends on the (unknown) radial distribution 
of the emitting clouds, this model does not constrain the 
FWHM/(Tiinc ratio. The thin disc model is able to explain 
the / values derived from C iv, with ranging between 5 
and 30 degrees. The /-FWHM relation given in figure [5] is 
accounted for by difi'erent inclination angles. On the other 
hand, H/3 data are suggestive of a wider range of ■& values, 
up to 90° in the thin disc picture. This mismatches with the 
Type-1 AGN model. Furthermore, objects with both UV 
and optical spectra exhibit difi'erent / values, for C iv and 
II/3 lines. Since / depends only on i9, different inclination 
angles for the orbits of the clouds emitting the two lines are 
needed. 

Geometrically thick disc model: 

In the case of a thick disc model C2 is non-negligible. Since 
H/3 is emitted in a larger region than C iv, the disc is thinner 
in the inner region (where C iv line i s emitted), and th icker 
outside (see the flared disc model in lCoUin et al. I l2006l . and 
the references therein). This model accounts for the differ- 
ences in the / factor of the two lines, and the non-correlation 
of the FWHMs of H/3 and Civ (figure 0. When the disc is 
seen almost face-on, the velocity component perpendicular 
to the disc plane would be larger than the projected ro- 
tational component, thus leading to FWHM values larger 
for H/3 than for Civ. This picture also explains why H/3 
FWHM/cTiinc ratios are close to the expected values for a 
thermal energy distribution at a given radius from the black 
hole. 
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APPENDIX A: 

Al Fell subtraction 

The reliability of our zero-order correction was checked as 
follows: 



(i) We observed I ZwOOl with the grism 7 setup. The H/3 
line was modeled and removed with a 2-gaussian fit. 

(ii) For each target, we applied the zero order correction 
for the Fe ll contamination, and we fitted the broad emission 
of H/3 with a 2-gaussian profile. In order to avoid the H/3 
narrow component, we excluded the line central region (~ 
1.5 times the spectral resolution) in the fitting procedure, 
when a narrow H/3 component was clearly observed. 

(iii) The spectrum of I ZwOOl has been convolved to 
a gaussian mimicking both the intrinsic and instrumental 
broadening of the Fe ll features. The line width was assumed 
to be the same as the one measured for H/3, consistently with 
the most of the literature, thus assuming that Fell and H/3 
emitting regions are nearly the same (but see iKuehn et al.l 
120081 ). 

(iv) The Fe ll template was then scaled in fiux and wave- 
length to match the observed features in the 4400-4650 and 
5100-5350 A ranges. The best fit was subtracted to the ob- 
served spectrum. 

(v) The narrow lines are modeled on the [O ml ° line, 

^ L 5007 A 

following iMcGill et al] l|2007l '): the [Olll] = fiux was as- 
sumed to be 1/3 of the [OmlgQQy^ line. Narrow H/3 and 
Hell fiuxes are left free in the fitting procedure. 

(vi) The resulting spectra only present the H/3 broad 
emission. We applied a 2-gaussian fit to the H/3 broad com- 
ponent, fitting the range 4750 - 4975 A. 

The template-subtracted FWHM and CTiino are compared 
with the one presented in the p aper in figure lAll u pper 
panels. As previously reported bv lMcGill et al.l |2003), the 
FWHMs obtained with and without the subtraction of the 
Fell template are in very good agreement: the average dif- 
ference in the two estimates is negligible (~ 0.03 dex), the 
standard deviation being ~ 11 per cent of the FWHM val- 
ues, comparable to the estimated error in the fitting proce- 
dure. The estimates of cruno are in agreement, too, but the 
dispersion is larger: the average difference is ~ 0.04 dex, and 
the residual standard deviation is ~ 20 per cent. 

A2 The fitting function 

Several authors (e.g., [McGill et al.l 12007') preferred the 
Gauss-Hermite series Dcr Marel & Franx 1993) when 

fitting the profile of broad emission lines. In this technique, 
the observed profile is fitted with a set of orthonormal poly- 
nomials (the Hermite ser i es) m ultiplied to a Gaussian curve. 
I Van Der Marel fc Fran^ (|l993h proved that, in most of the 
situations of astrophysical interest, the lines are well fitted 
with a series truncated at the fourth order. In this way, 
the set of independent parameters provided by the fit has a 
straightforward interpretation, in particular and h4 (the 
coefficient for the third and fourth order of the Hermite se- 
ries) are related to the line asymmetry and kurtosis. 

In order to check the dependence of our results on the 
adopted fitting function, we applied the Gauss-Hermite fit 
(extended to the fourth order of the series) to our data, after 
the Fe ll template subtraction described in appendix lAll The 
FWHM and crunc values obtained with the two profiles are 
compared in figure fAT] lower panels. The FWHM estimates 
are in good agreement with those adopted in the paper: 
the average difference (~ 0.05 dex) is negligible within the 
purposes of our work, with a residual standard deviation of 
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Figure Al. The comparison between various estimates of H/3 
FWHM (left) and uiine (right). In the x-axis, the width estimates 
obtained with our zero-order correction for Fe II emission and the 
2-gaussian fit for the broad emission of H/3. In the y-axis, in the 
upper panel the width estimates obtained with the 2-gaussian 
technique, after Fell template subtraction. In the lower panel, 
the width estimates based on the Gauss-Hermite fit proecdure. 
The dashed lines refer to the 1-to-l case. Typical error boxes are 
also shown. All FWHM estimates are consistent. As expected, the 
(Tiine values are more sensitive to the correction for Fe ii emission 
and to the fitted function. 



■~ 13 per cent. A systematic deviation of anno estimates is 
observed when comparing the two fitting technique, but the 
overall average difference (~ 0.08 dex) cannot account for 
the differences in the FWHM/criine ratio observed for CiV 
and H/3. 
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